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Abstract. We report on the search for B^ — > [i~ , B° — > /i fi~ decays and b — > sfi /i~ transitions 
in exclusive decays of B mesons using the CDF II detector at the Fermilab Tevatron Collider. Using 
2 fb -1 of Run II data we find upper limits on the branching fractions B(B^ — > /i + /i~) < 5.8 x 10~ 8 
and B(Bd — > < 1.8 x 10 -8 at 95% confidence level. The results for the branching fractions of 

the b -> sfi + fi~ transitions using 924 pb" 1 of Run II data are B(B + -> fi + fj,~K + ) = (0.60 ± 0.15 ± 
0.04) x lO" 6 , B(B° -> fi+fi-K* ) = (0.82 ±0.31 ±0.10) x HT 6 and B(B° -» fi+ n~ 4>)/B(B° s -> 
J/ip4>) < 2.61 x 10" 3 at 95% confidence level. 

PACS. 13.25.Hw Decays of bottom mesons - 14.40.Nd Bottom mesons 



1 Introduction 

The decay of a b quark into two muons, as well as in 
an s quark and two muons, requires a flavor-changing 
neutral current (FCNC) process which is highly sup- 
pressed in the standard model (SM) as they can only 
occur through higher order diagrams. New physics can 
significantly enhance the branching fractions of these 
decays. In these proceedings the current results of the 
CDF experiment are presented for the branching ratios 
of the rare decays -B( s( j) ~~ * an( l B ~ ¥ iJ, + fi~h, 

where B stands for B + , B®, or B s , and h stands for 
K + ,K*° or <j>. The K*° is reconstructed in the mode 
K*° — > K + n~ and the <j> is reconstructed as (f) — * 
K + K~ . A detailed description of the analyses can be 
found in Ref. [HE]. 



2 The CDF II Detector 

The CDF II detector is a cylindrical general-purpose 
particle detector built at one of the two collision points 
of the Tevatron pp collider which operates at a center- 
of-mass energy of = 1.96 TeV. Its inner track- 
ing system consists of a silicon microstrip detector 
surrounded by an open-cell wire drift chamber. The 
tracking systems are immersed in a 1.4 T magnetic 
field and measure the momentum of charged particles. 
The electromagnetic and hadronic sampling calorime- 
ters are located outside the solenoid. The outermost 
part of the CDF II detector is the muon detector sys- 
tem. Muons are detected by the planar drift cham- 
bers (CMU) and the central muon extension (CMX), 
which consists of conical sections of drift tubes. The 
CMU covers a pseudorapidity range up to \r)\ < 0.6, 
where r\ = — ln(tan | ) and 8 is the angle of the track 
with respect to the beamline, selecting muons with a 
Pt > 1.4 GeV/c. The CMX extends the coverage to 



a pseudorapidity range of 0.6 < |7y| < 1.0 for muons 
with p T > 2.0 GeV/c. 



3 Search Methodology 

The searches for the rare decays B® s ^ — > /.i + fi~ and 

B — > fi + fi~h use in both cases a similar approach. A 
data sample with an integrated luminosity of 2 fb _1 is 
used for B° s ^ — * n + search, respectively 924 pb -1 

for the B — * /i + fi~h search. In both cases events are 
selected by the dimuon trigger. For the B° s ^ — * n + 
analysis the data is futher divided into two classes. Ei- 
ther both muons are reconstructed in the CMU cham- 
bers, further called CMU-CMU, or one muon is recon- 
structed in the CMX chambers and the other in the 
CMU chambers, called CMU-CMX. 




Fig. 1. Graphical representation of the discriminating vari- 
ables. 
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3.1 Selection Optimization 



3.3 Background Estimation 



To optimize the data selection a signal event sample 
from MC simulation and a background sample from 
data sidebands is used. 

In case of the B® s « — ► ji + fi~ analysis a multivariate 

neural network (NN) enhances the signal and back- 
ground separation. It is based on the discriminating 
variables illustrated in Figure[TJ the proper decay length 
A = fljcT, the 3D opening angle AO between the dimuon 
momentum p w and the displacement vector from the 
primary to the dimuon vertex L, the pt of the lower 
momentum muon candidate and the B-candidate iso- 
lation / [TT] , 

In case of the B — > ^ /j,~h analysis a cut based ap- 
proach is used to optimize the selection for the figure 
of merit 5*/ ^ S + B, where S is the estimate of the ex- 
pected yield of the rare decays and B is the expected 
background. The discriminating variables are the de- 
cay length significance A/cta, the pointing angle a from 
the B meson candidate to the primary vertex and the 
isolation /. For the final selection candidates with a 
dimuon mass near the J/ip and the are rejected. 



3.2 Normalization Modes 

In order to obtain the branching ratios of the different 
rare decays, one normalizes to the B — > J/tjjh modes. 
The branching ratios for the decays B — > fi + [i~h are 
then given by 



B(B — > h) _ N^+p-h ej/^h 

B(B — > J/iph) Nj/^h e^+^-h 



B(J/^ 



n + n ) 

(1) 



where N,, 



is the observed number of B — ► h 



' ^ p. h 

decays, Nj/^h is the observed number of B — > J/iph 
decays, while ej/^h an d c^+^-ft are the selection effi- 
ciencies of B — > J/4>h and B — > fi + [i~h respectively. 
The ratio of efficiencies is about 70 to 85 % [JJ. 
In case of the rare decays B® s d ^ — > the upper 

limit on the branching fraction can be expressed as 
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In case of the B® d 



l l+ 1 1 



consists of contributions from B® d 



To estimate the background contributions several dif- 
ferent sources are considered. 

In case of the rare decays B — > n + n~h these sources 
are charmless B decays into charged hadrons, reflec- 
tions between the three rare decay modes and combi- 
natorial background. The background originating from 
charmless B decays is calculated from a simulation 
of these decays convoluted with the misidentification 
rates of muon detectors obtained from Z?*-tagged D° — > 
K~tt + decays. The combinatorial background is es- 
timated from the high mass sidebands of the B sig- 
nal and extrapolated under the signal region using 
the background shape from the data distribution with 
poor vertex quality. 

decays the background 
h + h~, where 

h = tt^ or K^, and combinatoric background. The 
expected background from B® d — > h + h~ is calculated 
from equation[2]by replacing B® d — > with B® d — > 

hrh~ and including two additional efficiency terms to 
account for the muon misidentification rates, measured 
in D° — ► Kir data, and the fraction of misidentificd 
B® d — > h + h~ events falling in the signal windows. The 
branching ratio for the various B® d — ► h + fi~ modes 
are taken from Ref. [3] . The combinatoric background 
is estimated by extrapolating the number of events in 
the sideband regions passing a given cut to the signal 
region using a linear fit. The total background is then 
formed by summing up the combinatoric background 
and the contributions from B° s d — > h + h~ decays. 



4 Results 



4.1 B 



H + Il- 



ia case of the rare decays B — > ii + /i~h we observe an 
excess in the signal region in all three decay modes. 
The invariant mass distributions are shown in Figures 
[2 [3] and |H The significance of the excess is determined 
by calculating the probability for the background to 
fluctuate into the number of observed events. The re- 
sults for the number of observed events, expected back- 
ground events, the significance and the absolute and 
relative branching ratios are listed in table [JJ The 
branching ratios are calculated using equation [TJ We 
find 



where N^ 5 % is the upper limit on the number of B s d — > 

H + l-i~ decays at the 95% C.L. determined from the 
comparison between expected and observed background 
events, Ng+ is the number of reconstructed B + — > 
J/tpK + candidates. The parameters a, e basc and e NN 
are the trigger acceptances and the efficiencies of the 
initial, respectively the NN, requirements. -4^- denotes 
the ratio between the probabilities that a b quark pro- 
duced in a ppbar collisions hadronizes into a B + meson 
and a B s or Bd meson [2]. 



B(B+ -> [i+ii-K + ) 
B(B+ -> J/tpK+) 
B(B° d -» fi+fi-K* ) 
B(B° d J/iiK* ) 
B{B a s -> 
B(B° -> J/U) 



= (0.59 ±0.15 ±0.03) x 10~ 3 
= (0.62 ±0.23 ±0.07) x 10~ 3 
= (1.24 ±0.60 ±0.15) x 10~ 3 



Using the world average branching ratios of the nor- 
malization modes [3], we find 

B{B + -> fi + fi-K+) = (0.60 ± 0.15 ± 0.04) x 10~ 6 
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B{B° d -> fi+fi-R* ) = (0.82 ± 0.31 ± 0.10) x 10~ 6 
B(B° S -> fi + fj,-(j>) = (1.16 ± 0.56 ± 0.42) x 10~ 6 

Since the excess in B® — > ji + /i~ <f> is not significant, we 
calculate a limit on its relative branching ratio using 
Bayesian integration assuming a flat prior, and find at 
95(90)% C.L. 
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Fig. 2. Inv. mass spectrum for the decay 
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Fig. 3. Inv. mass spectrum for the decay 
(i + fj,~ K*° } 
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Fig. 4. Inv. mass spectrum for the decay B s 
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In case of the rare decays B° s d — ► H + H~ we use dif- 
ferent neural network bins and mass bins for the com- 
putation of the limits on the branching ratios. Table 
[5] shows the number of expected and observed events 
of the two trigger scenarios CMU-CMU and CMU- 
CMX for different cuts on the network output NN and 
gives the result on the branching ratio for the com- 
bination of both scenarios with a neural network cut 
of NN > 0.80. In Figures [5] and [5] the invariant mass 
distribution vs. the neural network value is displayed. 
Using equation [2] we obtain the 90(95)% C.L. limits 
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Fig. 5. The invariant mass distribution versus the NN 
output for the CMU-CMX channel. 



CDF II Preliminary 



(2 fbjj 




Fig. 6. The invariant mass distribution versus the NN 

1 The solid line is a graphical representation of the ex- output for the CMU-CMU channel, 
tracted yield, not a fit to the data. 
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Table 1. Results of the B — > (J, fi h analysis. 



Mode 


B+ -> [1+H-K+ 


B d -» /i /i K 


B° — > 


N s 


90 


35 


11 


N BG 


45.3 ±5.8 


16.5 ±3.6 


3.5 ± 1.5 


Gaussian significance (a) 


4.5 


2.9 


2.4 


Rel B ± stat ± sys x 1CT 3 


0.59 ±0.15 ±0.03 


0.62 ±0.23 ±0.07 


1.24 ±0.60 ±0.15 


Abs B ± stat ± sys x 1CT 6 


0.60 ±0.15 ±0.04 


0.82 ±0.31 ±0.10 


1.16 ±0.56 ±0.42 


Rel B 95% CL limit x 1(T 3 






2.61 


Rel B 90% CL limit x 1(T 3 






2.30 



Table 2. Expected number of background (exp.) and number of observed (obs.) events in the CMU-CMU and CMU- 
CMX channels for the B° mass window (5.310-5.430 GeV/c 2 ) and the B° d mass window (5.219-5.339 GeV/c 2 ). 





Mode 


B° d -» /x+m" 


B° s 


CMU-CMU 


N exp., 0.800 <NN< 0.950 


25.5 ±0.7 


23.5 ±0.7 


CMU-CMU 


N obs., 0.800 <NN< 0.950 


32 


18 


CMU-CMU 


N exp., 0.950 <NN< 0.995 


8.5 ±0.5 


7.7 ±0.5 


CMU-CMU 


N obs., 0.950 <NN< 0.995 


7 


10 


CMU-CMU 


N exp., 0.995 <NN 


2.4 ±0.2 


2.1 ± 0.2 


CMU-CMU 


N obs., 0.995 <NN 


5 


2 


CMU-CMX 


N exp., 0.800 <NN< 0.950 


27.8 ±0.9 


26 ± 0.7 


CMU-CMX 


N obs., 0.800 <NN< 0.950 


28 


26 


CMU-CMX 


N exp., 0.950 <NN< 0.995 


10.8 ±0.5 


10.3 ±0.5 


CMU-CMX 


N obs., 0.950 <NN< 0.995 


6 


11 


CMU-CMX 


N exp., 0.995 <NN 


1.6 ±0.2 


1.6 ± 0.2 


CMU-CMX 


N obs., 0.995 <NN 


1 


1 


combined 


Abs B 95% CL limit xlO -8 


1.8 


5.8 


combined 


Abs B 90% CL limit xlO" 8 


1.5 


4.7 



5 Conclusion 



We present the latest results on the measurement of 
the branching ratios of the rare decays B — > [i + fi~h 
and give an upper limit on the branching ratio of the 
decays B° s d — > The results for the branching 

ratios for the decay modes B + — > [i + ^~ K + and £?° — > 
fi + fi~K*° are in good agreement with the results of 
the B factory experiments BABAR and BELLE [9lfT0]. 
The limit on the branching fraction B(B® — * 
S(i?° — * J/ip(f>) is the most stringent to date. The new 
results for the upper limit for the branching ratio of 
the decays B° s d — * ^ + ^~ are currently the world's best 
limits and can be used to reduce the allowed parameter 
space of a broad spectrum of SUSY models [5,6j3H]. 



6 Acknowledgments 

The author would like to thank the members of the 
CDF Collaboration who performed the analyses. 



References 

1. CDF Collaboration, CDF Public Note 8543, (2006) 

2. CDF Collaboration, CDF Public Note 8956, (2007) 

3. A. Abulencia et a/., Phys. Rev. Lett. 97, 211802 
(2006) 

4. W.M. Yao et al, J. Phys. G 33, 1 (2006) 

5. R. Dermisek et al, J. High Energy Phys. 09, 029 
(2005) 

6. R. Ruiz de Austri et al., J. High Energy Phys. 0605, 
002 (2006) 

7. S. Baek et al, J. High Energy Phys. 0506, 017 (2005) 

8. R. Arnowitt et al, Phys. Lett B 538, 121 (2002) 

9. B. Aubert et al. [BABAR Collaborat ion], Ph ys. Rev. 
D 73, 092001 (2006), [arXlv:hep-ex/0604007| 

10. K. Abe et al. [BELLE Collaboration], 
|arXiv:hep-ex/0410006| 

11. The B-candidate isolation is defined as I — 
IPt M Pt ± |p T |)> where the sum is over all tracks 
with yj An 2 + A<j) 2 < 1; Acj> and Arj are the azimuthal 
angle and pseudorapidity of track i with respect to 



